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ABSTRACT 

I describe the principle of using the 1.6/iin H~ spectral feature as a photometric 
redshift indicator and demonstrate that the technique holds promise by successfully 
recovering the redshifts of a small sample of 2=0-1 galaxies using only their infrared 
(JHKL) photometry. I then consider the applicability of the technique to the 3.6-8/im 
SIRTF filter set and investigate the systematic errors that could arise in photometric 
redshifts from random photometric errors or from a mismatch between target galaxies 
and fitting templates in metallicity, star formation history, and amount of interstellar 
dust. It appears that SIRTF near-IR data alone should be sufficient to estimate redshift 
of most galaxies that are at 2>1.5 and are dominated by stellar populations older 
than >20Myr. Galaxies whose photometric fits indicate them to be at lower redshifts, 
2fit<1.5, or dominated by very young stellar populations, agejj t <20Myr, suffer from 
severe degeneracies in photometric redshift, and a reliable photometric determination 
of their redshifts will have to include either IR observations at shorter wavelengths 
(H and K) or optical data. Overall, it appears that — with care and caveats — the 
1.6/iin bump can provide a powerful way of estimating redshifts of distant galaxies in 
deep infrared imaging surveys that will soon be provided by SIRTF, and, eventually, by 
NGST. 

Subject headings: galaxies: high-redshift — galaxies: distances and redshifts — tech- 
niques: photometric 



1. INTRODUCTION 

The study of galaxy formation and evolution will gain an enormous boost from space-borne 
observations in the infrared (IR), first with the advent of the Space Infrared Telescope Facility 
(SIRTF) and then, at the end of this decade, with the appearance of the Next Generation Space 
Telescope (NGST). SIRTF, scheduled for launch in early 2003, will allow us to efficiently observe 



- 2 - 



distant galaxies at mid- and far-infrared wavelengths. In particular, SIRTF's Infrared Array Camera 
(IRAC) will operate at 3.6/xm, 4.5/xm, 5.8/um, and 8/xm, providing a window on the rest-frame near- 
IR properties of galaxies at z>l. A number of deep SIRTF imaging surveys that will detect many 
thousands of galaxies to high redshift are already scheduled (e.g., GOODS 1 , SWIRE 2 ), and it will 
be enormously important to efficiently estimate redshifts of galaxies in these samples. However, 
large-scale spectroscopic follow-up of these IR-selected surveys will be very time-consuming and 
only small sub-samples of objects are likely to be targeted. Instead, spectroscopic redshifts are 
likely to be complemented by the less precise but much more efficient photometric redshifts. 

Photometric redshifts enjoy great popularity as a way to estimate the redshifts of faint galaxies 
(e.g., Connolly et al. 1995; Gwyn & Hartwick 1996; Sawicki, Lin, & Yee 1997; Fernandez-Soto, 
Lanzetta, & Yahil 1999; Rudnick et al. 2001). Traditionally, photometric redshifts rely on the 
photometric signatures of the Lyman and 4000A breaks to estimate redshifts from optical and 
near-IR data. However, the availability of mid-IR data from SIRTF (and, later, NGST) surveys 
opens up the possibility of using a hitherto untried photometric redshift technique that could 
obviate the need for optical imaging of SIRTF and NGST survey fields, or provide independent 
redshift estimates where such imaging is already present. Specifically, the spectral "bump" at 
1.6/um, caused by the minimum in the opacity of the H~ ion present in the atmospheres of cool 
stars, is a nearly ubiquitous feature of almost all stellar populations and the photometric signature 
of this bump can be expected to provide a means of estimating redshifts to distant galaxies. 

The importance of the imprint of the 1.6^m bump on SIRTF IRAC colors has been appreciated 
for a long time (e.g., Wright, Eisenhardt, & Fazio 1994), and indeed the filter transmission curves 
of the IRAC filters have been designed partially with photometric redshifts in mind by modeling 
the expected colors of high-z galaxies (Simpson & Eisenhardt 1999). The purpose of the present 
paper is to demonstrate the potential of the 1.6//m bump as a photometric redshift indicator and 
to explore some of the possible limitations of the technique at high redshift. Section 2 illustrates, 
with real data at 0<z<l, that the 1.6/xm bump is an effective photometric redshift indicator and 
describes the source of the photometric redshift signal. Section 3 discusses the applicability of the 
technique to SIRTF IRAC data and explores some of the limitations that could be present at high 
redshift due to effects such as metallicity, dust, and star formation history. The main conclusions 
of the paper are summarized and discussed in Section 4. 



see http: / /www. stsci.edu /science /goods / 
2 see http://www.ipac.caltech.edu/SWIRE/ 
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2. PRINCIPLE AND PROOF OF CONCEPT 

2.1. The 1.6^m Bump 

The opacity of the H~ ion has a minimum at ~1.6/xm (e.g., John 1988), which imprints itself as 
a maximum, or "bump," on the spectral energy distributions (SEDs) of cool stars. Consequently, 
the SEDs of composite stellar populations that contain significant numbers of such stars can be 
expected to possess the same spectral imprint. This assertion is illustrated in Fig. 1, which shows 
the spectral energy distributions of model stellar populations from the 1996 version of the GISSEL 
spectral synthesis code of Bruzual & Chariot (1993). The 1.6^m bump is easily discernible in 
all but the extremely young (~1 Myr) stellar populations, where it is swamped by the essentially 
power-law spectra of very hot, massive young stars. Given the strength and near universality of 
the 1.6^m bump, it is reasonable to expect that it will imprint a significant signal on the colors 
of distant galaxies, and hence, as was pointed out by Simpson & Eisenhardt (1999), that it might 
make a good photometric redshift indicator. 

2.2. The 1.6^m Bump in the Hubble Deep Field 

While the 1.6/xm bump has been suggested as a photometric redshift indicator before, its 
feasibility remains to be tested with real data. The JHKL filter set straddles the 1.6/xm bump 
at z«0-l, and so JHKL photometry of Hubble Deep Field (HDF; Williams et al. 1996) galaxies 
with spectroscopic redshifts will be used to test if the 1.6/xm bump can provide a strong enough 
photometric imprint to constrain redshifts. Note that the JHKL filter set provides a reasonable 
surrogate of the IRAC quartet of bandpasses, albeit shifted by a factor of ~2.5-3 in observed 
wavelength: while filters of the JHKL set cover the rest-frame 1.6/«n bump for galaxies at 
the corresponding redshift range for the IRAC filters is zpsI.8-4. Consequently, JHKL photometric 
redshifts of z<l galaxies provide a useful demonstration of how well SIRTF may do at these higher 
redshifts. 

Nine galaxies were detected at L-band (3.2/xm) in the deep Keck imaging of the HDF by 
Hogg et al. (2000). All nine have redshifts at 0<z<l, as measured by Cohen et al. (1996, 1999) 
and Lowenthal et al. (1997). The L-band photometry of these objects was adapted from Hogg 
et al. (2000), while their JHK magnitudes were measured from the public Kitt Peak infrared 
images of the HDF (Dickinson, 1998; Dickinson et al. 2002, in prep.). The JHK magnitudes were 
measured in 2"apertures using the SExtractor package (Bertin & Arnouts, 1996), while the L-band 
magnitudes of Hogg et al. (which they measured in 2"apertures but then extrapolated to 6" ones) 
were transformed back to 2" apertures by scaling (in reverse of the procedure used by Hogg et al.) 
the L-band 6" magnitudes by the ratio of LT-band fluxes measured in 2" and 6" apertures. The 
resulting 2"-aperture magnitudes are listed in Table 1. 

Photometric redshifts of the nine galaxies were estimated by comparing their observed broad- 
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band spectral energy distributions (as listed in Table 1) with model templates generated from the 
1996 version of the GISSEL spectral synthesis package (Bruzual & Chariot 1993). Solar-metallicity 
GISSEL spectra of constant star formation rate (SFR), and with ages spanning the full available 
range of models (0-20 Gyr), were used as the starting point in the construction of the templates. 
These model spectra were redshifted onto a grid of redshifts spanning z=0-1.5 in steps of 5z=0.05, 
and were then integrated through the filter transmission curves (FTCs) to produce the template 
broadband fluxes 



in the four filters i=(J, H, K, L) . The resulting grid contains 6851 templates, each of four broadband 
fluxes, and spans 0<z<1.5 linearly and 0<age<20 Gyr quasi-logarithmically. 

For each observed object, the observed fluxes were compared to each template in the template 
grid by computing the goodness-of-fit measure 



where f bs{i) and a(i) are the observed flux and its uncertainty in the i th filter, ft p t(i) is the flux 
of the template in that filter, and s is the scaling between the observed and template fluxes that 
can be computed analytically using 



For each object, the best-fitting template — and, hence, the most likely redshift — is then identified 
as the one that produces the smallest value of % 2 . Note that variants of this procedure are a 
common way of computing photometric redshifts using model templates (e.g., Gwyn & Hartwick 
1996; Sawicki et al. 1997, Arnouts et al. 1999), with the key difference being that the data and 
templates used here do not span the traditional 4000A or Lyman breaks, but, instead, cover the 
1.6/um bump. 

The fidelity of the photometric redshifts obtained using solely JHKL photometry is remarkably 
good, as is shown in Fig. 2 which compares the photometric and spectroscopic redshifts of the nine 
HDF galaxies (Table 1 lists the z p hot values object by object.). The scatter between photometric 
and spectroscopic redshifts is only o"a z =0.09, which is similar to the scatter typically obtained over 
the same redshift range in the HDF using the Balmer-break sensitive optical photometric redshift 
techniques (e.g., Sawicki et al. 1997; Fernandez-Soto et al. 1999). Although the comparison between 
photometric and spectroscopic redshifts is based on the very small sample of nine galaxies, the tight 
correlation between z p hot and z spec demonstrates, for the first time, that the 1.6^m bump is a viable 
indicator of redshift. 




(1) 
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2.3. The Source of Photometric Redshift Signal 

To understand how the l.Qfim bump produces a photometric redshift signature, it is instructive 
to consider how the infrared colors of galaxies depend on redshift. The JHKL filter set used in 
§2.2 can be represented as six constituent colors; Fig. 3 shows these six colors for GISSEL models 
of constantly star-forming stellar populations of seven different ages. The solid lines represent 
colors of stellar populations ranging in age from 20 Myr to 20 Gyr, while the dotted line is for the 
extremely young 1 Myr-old population. 

For the extremely young, 1 Myr-old model, all six colors are virtually constant with redshift: 
the spectra of the very young, hot stars that dominate this ultra-young stellar population arc 
essentially featureless Rayleigh-Jeans power laws (see Fig. 1) and because redshifting a power 
law does not affect its spectral shape, the colors of extremely young stellar populations remain 
degenerate with redshift. Consequently, the redshift of a galaxy dominated by an extremely young 
stellar population cannot be determined photometrically. Fortunately, stellar populations older 
than a few Myr develop the 1.6/um bump (see Fig. 1), and it is the passage of this bump through 
the JHKL filter set that imprints itself on the colors of galaxies in Fig. 3. For the remainder of 
this section, extremely young (~1 Myr-old) stellar populations will be ignored and the analysis will 
focus on the more prosaic models with age>10 Myr. 

Clearly, colors which are degenerate with redshift cannot provide redshift information: for 
example, the J—H colors of all model ages are essentially flat with redshift and so J—H contains 
little information about redshift. In contrast, where a color changes steeply with redshift, such as 
K—L over 0.4<z<0.8, it provides strong redshift signal. However, to be useful, the color must not 
only change steeply, but must also be unique with redshift: thus, although the J—L color changes 
steeply with redshift at J— L=— 0.5, it does not constrain redshift well since models of different age 
produce that same color at different redshifts. Therefore, the requirements for a color to be a good 
photometric redshift indicator are that (1) it change steeply with redshift and (2) it be unique for 
a given redshift. 

Fig. 3 shows that J—H is a poor redshift indicator over all redshifts z<1.5; H—K is fine only 
at z~0; J—K would be fine at low redshift were it not for the fact that in this color z«0 objects 
are degenerate with young stellar populations at z>l; J—L is fine at z&0, H—L at 0<z<0.4, and 
K—L at 0.4<z<0.8. Thus, the success of the JHKL filter set at determining the redshifts of the 
nine HDF galaxies is primarily due to the H—L color below z~0.4 and K—L at higher z. 

The two most useful colors, H—L at z<0A and K—L at z>0A, have an important feature 
in common: they both contain the longest-wavelength L filter. The role of this filter may be 
appreciated in Fig. 1: At all redshifts z<l, the L filter is redward of the 1.6/um bump — i.e., at 
wavelengths at which the spectra of virtually all stellar populations are nearly identical. Thus, the 
L filter provides a universal "anchor," which is independent of galaxy type. The role that the other 
two filters (H and K) play now becomes apparent: while L acts as a redshift-invariant anchor, the 
1.6/rai bump moves with increasing redshift first through J and then through K. The passage of 



-6- 



the 1.6/xm bump through the J and K bandpasses depresses the fluxes in these filters with respect 
to the invariant L-band anchor, thus resulting in the steep yet age-invariant evolution of the H—L 
and K—L colors over 0<z<0.4 and 0.4<z<0.8, respectively. 

As is shown in Fig. 1, spectra of galaxies are virtually identical redward of rest-frame 1.5^m and 
so the H—L and K—L colors of galaxies are independent of age over the, respectively, 0<z<0.4 and 
0.4<z<0.8 redshift ranges. At wavelengths below rest-frame 1.5//m, galaxy colors become strongly 
dependent on age and, consequently, at z>0A, the H—L color develops an age-redshift degeneracy. 
However, the K—L color takes over at z>0.4 as a redshift indicator, and the H—L color now can 
assume a different role — that of providing information about the evolutionary state of the galaxy's 
stellar population! 

The above discussion illustrates the general principle of selecting filters for measuring redshifts 
using the 1.6//m bump: One needs an "anchor" filter which will at all redshifts of interest remain 
in the age-invariant part of the spectrum redward of the 1.6/xm bump. Furthermore, one needs 
a filter whose bandpass corresponds to the wavelengths through which the 1.6//m bump will be 
passing in the redshift range of interest. A third filter can then be deployed in the age-sensitive 
region blueward of the 1.6/zm bump to gain information about the state (age) of the galaxy's stellar 
population. Finally, to cover an extensive redshift range, a set of adjacent filters is necessary so that 
the filters can take over from each other at different redshifts. In practice, of course, fluxes from 
all available filters contribute when the photometric redshift is determined via x 2 minimization of 
broadband spectral energy distributions; nevertheless, the bulk of the photometric redshift signal 
at a given redshift z comes from the combination of a long-wavelength anchor bandpass and a 
bump-containing filter at \futer~(,l + z)1.6//m. 



3. THE l^m BUMP AND SIRTF 

Armed with the understanding of how the 1.6//m bump interacts with broadband filters, let 
us consider how the SIRTF IRAC 3.6/xm, 4.5/xm, 5.8^m, 8/xm bandpasses can be used to estimate 
redshifts. The 8/im bandpass is the reddest of the IRAC filters and so it naturally serves as the 
long-wavelength anchor. In combination with the 8/xm anchor, the bluer 3.6/um, 4.5/xm, and 5.8fim 
bandpasses then provide the bulk of the photometric redshift signal over different redshift ranges. 
This principle is illustrated in Fig. 4, which shows the IRAC colors of model stellar populations 
(see also Fig.9 of Simpson & Eisenhardt 1999). The ni3.6— m§ color provides a good photometric 
redshift indicator over 1.3<z<2, while becoming degenerate to galaxy age at higher redshifts; the 
11L4.5— nig color works at 2<z<2.7, while ms.g— mg is good at 2.7<z<3. 

Note that at low redshift (z<l-2, depending on the central wavelength of the blue filter) 
IRAC filters are not probing the region around the 1.6/um bump. Instead, they are sensitive to 
the presence of the CO absorption band at 2.4/xm. The interplay of the IRAC filters with the CO 
band results in a series of "wiggles" in the colors at z<l-2. These wiggles make the IRAC colors 
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degenerate at low redshift and so restrict the redshift range over which IRAC photometry can be 
used to measure redshift. 

Let us now consider in more detail what biases and uncertainties can be expected when trying 
to determine redshifts using SIRTF IRAC photometry. Photometric redshift measurements can 
be biased, sometimes catastrophically, by random errors in photometry, as well as, in a more 
systematic way, by incorrect assumptions about the target galaxy's spectral energy distribution. 
Section 3.1 will examine the biases that can be introduced by simple Gaussian scatter in photometric 
measurements, while §3.2 will examine the problems that may arise from incorrect assumptions 
about the metallicities, star formation histories, and dust content of high-z galaxies. 

3.1. Random Photometric Errors 

In practice, photometric redshifts of SIRTF-selected galaxies will be computed using all four 
IRAC bandpasses (plus possibly additional observations at other wavelengths), using a method 
similar to the template-fitting technique that was used to determine photometric redshifts of HDF 
z<l galaxies in §2.2. To understand how random Gaussian photometric errors can affect the 
measurement of photometric redshifts, it is instructive to consider how such photometric errors 
propagate through eq.(2). 

To simulate the IRAC colors of high-z galaxies, a grid of models was computed based on the 
constant star formation GISSEL spectra covering the full range of available ages 0<age<20Gyr 
and spanning 0<z<5 in steps of fc=0.05. Since high-z galaxies can be expected to have sub-solar 
metallicities, Z=0.2 Z & GISSEL SEDs were used. These SEDs were integrated through the four 
SIRTF IRAC filter transmission curves (see eq.[l]) to produce a grid of predicted fluxes. Sample 
"observed" galaxies were taken from this model grid and were fit via eq.(2) with the full grid of 
models as the template grid. Each "observed" galaxy was assigned photometric errors of cr(i)=(0.05, 

0. 05, 0.75, 0.1) in the four IRAC filters, z=(3.6^m, 4.5^m, 5.8/um, 8/mi), which correspond to the 
errors expected for a z=3 L* galaxy in the GOODS fields. 

Fig. 5 shows the fit results in redshift-age space for ten "observed galaxies" — at two input, 

1. e. "true," ages (agej„=50Myr in the left-hand panels and lGyr in the right-hand panels) and five 
redshifts, Zj n =0.5, 1, 2, 3, and 4 (top to bottom). In each panel, the true redshift-age location of 
the "observed" galaxy is shown as a filled circle; since the "observed" galaxies are taken directly 
from the fitting grid, this is also, by construction, the location of the best-fit template identified by 
the fitting process. Regions of fit confidence are identified by selecting those templates for which 
X 2 values depart by less than a certain amount from the best-fitting Xmin'i shown are regions of 
95% confidence for "observed" data with random Gaussian photometric errors of cr(i)=(0.05, 0.05, 
0.075, 0.1) in the IRAC filters i=(3.6fim, 4.5/im, 5.8/um, 8/mi). It is the shapes of these confidence 
regions that illustrate the biases and scatter that will exist in photometric redshifts derived from 
SIRTF IRAC data. 
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The inspection of Fig. 5 leads to a number of observations: 

1. The girth of the contours in the top two rows of Fig. 5 is indicative of the expected scatter in 
photometric redshifts at z<1.5 and illustrates that the determination of photometric redshifts 
with IRAC data alone is unfeasible at z<1.5. This large scatter within 0<z<1.5 results 
from the degeneracies in IRAC colors seen at low z in Fig. 4. Thus, IRAC colors alone are 
insufficient for determining redshifts of z<\.5 galaxies with any precision. However, following 
the principles of filter selection that were detailed in §2.3, it should be possible to determine 
accurate photometric redshifts of z<1.5 galaxies in SIRTF surveys by augmenting their IRAC 
fluxes with H- and .ftT-band data. Of course, if optical data are available — as, for example, 
will be the case for the SIRTF GOODS survey — it will be possible to also use "traditional" 
(4000A-break-dependent) photometric redshifts to estimate reliable distances to z<l SIRTF- 
selected objects. 

2. In addition to the large scatter within 0<z<1.5, the contours in the top two rows of Fig. 5 show 
that z<1.5 galaxies can be easily mistaken for very young (age <10Myr) stellar populations at 
nearly any redshift. Conversely, extremely young stellar populations at high redshift (z>1.5) 
can easily be misidentified as z<1.5 galaxies. This degeneracy is a result of the similarity 
between the colors of low- z objects and those of very young stellar populations at any redshift 
(Fig. 4). Because of this degeneracy, SIRTF data alone cannot be used as a reliable way to say 
that a galaxy is at z<1.5; neither can they be used to unambiguously identify high-z galaxies 
that are dominated by very young stellar populations. Such degeneracies could possibly be 
broken with the addition of UV-optical data that enable the use of the Lyman break technique 
(e.g. Steidel et al. 1996) or of "traditional" photometric redshifts (e.g. Sawicki et al. 1997). 

3. Photometric redshift of galaxies whose fits show them to be both at z>1.5 and older than 
~20 Myr can be relied on: there exists no confusion in Fig 5 between z<1.5 galaxies and 
older higher- z objects, at least under the assumption (in force here) that all the errors in 
photometric redshifts arise from random errors in photometry. Rejecting galaxies whose 
photometric redshifts claim them to be either at z<1.5, or younger than ~20 Myr, leaves a 
reliable core of objects which should be free of catastrophic errors in photometric redshift. 

4. The expected scatter in photometric redshifts of objects with z>1.5 and ages>20Myr is 
smallest at z^2 and increases to higher redshifts. This increase results because IRAC colors 
become progressively flatter and more degenerate with increasing redshift (see Fig. 4) as the 
1.6//m bump encroaches into the reddest (anchor) bandpass. The addition of a redder filter 
would provide an anchor bandpass that would improve the accuracy of photometric redshifts 
at z~3 and beyond. Fortunately, with IRAC data alone, the scatter is small at z~2, which 
is precisely the epoch where spectroscopic redshifts are hardest to secure. The epoch around 
zKi2 is also the era where the cosmic star formation appears to plateau (e.g., Sawicki, Lin, & 
Yee 1997; Steidel et al. 1999) after its strong rise with lookback time from z=0 to z«l (Lilly 
et al. 1996); the ability of SIRTF data to reliably select z^2 galaxies makes SIRTF surveys 
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well-poised to fill the gap in our picture of the galaxy population at this crucial transitional 
epoch. 

5. IRAC data alone are insufficient to provide strong constraints on galaxy ages. The ages of 
intrinsically very young stellar populations (age<10Myr) could in principle be constrained 
well if one could break the strong degeneracy between such objects and those at z<1.5 by 
the addition of bluer filters, or — ideally — spectroscopic redshifts. The ages of intrinsically 
older galaxies (age>50 Myr) are more difficult to constrain because IRAC colors at high z do 
not evolve quickly with time. However, while the scatter in photometric redshifts increases 
with redshift past z«2, the accuracy of age determination improves because of the divergence 
of (primarily) UI3.6— m§ colors at higher redshifts; the addition of bluer IR filters should allow 
ages to be constrain even more precisely. 

Thus, it appears that SIRTF IRAC colors alone may be sufficient to measure photometric 
redshifts of high-z galaxies, provided that care is taken to reject objects whose fits claim them to be 
either at z<1.5 or younger than ~20 Myr (these redshift and age cut-offs can be adjusted depending 
on the size of the photometric uncertainties and on the desired balance between completeness 
and contamination of the sample). At z<1.5 the accuracy of photometric redshifts would greatly 
improve with the addition of H and K fluxes, while at higher redshifts (z>2.5) low- 2 contaminants 
can potentially be separated from young high-z galaxies by the addition of UV-optical data that 
are sensitive to the presence of the Lyman break. Overall, random photometric errors of reasonable 
size appear to allow photometric redshifts to be determined using SIRTF data over selected regions 
of parameter space. 

3.2. Systematic Errors 

So far only the effects of random photometric errors have been considered. However, given 
our incomplete state of knowledge about the nature of high-z galaxies, there exist many potential 
sources of systematic error. If the spectral energy distributions of high-z galaxies are different from 
what we assume, our photometric estimates of these galaxies' redshifts may be biased. Metallicity 
effects, deviations from the expected makeup of the stellar populations, and the presence of dust can 
all modify the shape of the spectrum around 1.6//m, and thus can systematically bias photometric 
redshifts. This section considers what happens when the dust-free, constantly star-forming, Z=0.2 
Z & template grid of §3.1 is used to fit "observed" galaxies of discrepant metallicity, star formation 
history, or dust content. 
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3.2.1. Dust 

High-z galaxies can be shrouded by substantial amounts of interstellar dust. For example, 
Sawicki & Yee (1998) found that some Lyman Break Galaxies in the HDF may be reddened by as 
much as E(B — V)^0A, or Ay^l (see also Papovich, Dickinson, & Ferguson 2001; Shapley et al. 
2001). Because interstellar dust tends to preferentially suppress flux in the blue, it can be expected 
to both shift the apparent wavelength of the 1.6^m bump and to make the stellar population appear 
older. 

To study the effect that dust may have on the accuracy of photometric redshifts, templates 
constructed from constantly star-forming, Z=0.2 Zq, dust-free GISSEL models (i.e., identical to 
those in §3.1) were used to fit "observed" galaxies that had the same metallicity and star formation 
history as these fitting templates, but were reddened with E(B — F)=0.3 following the Calzetti 
(1997) starburst dust curve. Fit results are shown in Fig. 6; as with Fig. 5, input galaxies of 
agej n =50 Myr and lGyr, and located at Zj„=0.5, 1, 2, 3, and 4, were used. The asterisks in Fig 6 
indicate the locations of these input parameters, while — as in Fig. 5 — filled circles and contours 
show the best-fit values and regions of 95% confidence for <r(i)=(0.05, 0.05, 0.075, 0.1) Gaussian 
photometric errors. 

The dust-free template grid recovers faithfully (certainly to well within the 95% error bars 
in Zfit) the redshifts of the dusty "observed" galaxies. At z>3, the recovered ages are noticeably 
older than the input ones — presumably because dust preferentially suppresses flux blueward of 
the 1.6^m bump, thereby mimicking the redder SEDs of older stellar populations. Nevertheless, 
the ability to determine redshifts from IRAC photometry is not compromised by a mismatch in 
dust content between the target galaxy and the fitting template grid. Provided the earlier caveats 
against accepting fits with z^<1.5 or agejjt<20 Gyr are heeded, IRAC colors should allow faithful 
measurements of redshifts even in the presence of substantial amounts of dust. 

3.2.2. Star Formation History 

The star formation histories of high-z galaxies are at present not well understood (see, e.g., 
Sawicki & Yee 1998; Papovich et al. 2001; Shapley et al. 2001). Two extreme illustrative scenarios 
are (1) an instantaneous burst of star formation and (2) a constantly star-forming system. The 
effect of star formation history on the fidelity of photometric redshifts was studied by using the 
standard grid of constantly star-forming, Z=0.2 Zq, dust-free templates to fit "observed" galaxies 
that have the same metallicity and dust content as the grid, but which underwent only a single, 
instantaneous burst of star formation that occurred either 50Myr or lGyr (left and right-hand 
panels in Fig. 7, respectively) before being "observed." 

As Fig. 7 shows, despite a mismatch in star formation history between "observed" objects 
and fitting templates, redshifts of Zj„>1.5 galaxies are recovered faithfully. At Zj n >3, the ages of 
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galaxies tend to be overestimated — a result of the deficit of young, blue stars in the instantaneous- 
burst "observed" galaxies that makes them seem older than they are when fit with constant SFR 
templates. Nevertheless, IRAC-based photometric redshifts can be trusted provided that, as before, 
one avoids objects with z^<1.5 or age/jt^20Gyr. 

3.2.3. Metallicity 

The creation of the ion relies on the presence of a supply of free electrons, and so the 
abundance of metals in stellar atmospheres may affect the strength of the 1.6/xm bump. A metal- 
poor object can be expected to have a less pronounced 1.6^m bump and a modified shape of the 
SED either side of the bump, potentially affecting both the recovered redshift and age of the target 
galaxy. Metallicities of high-z galaxies are as yet poorly constrained, ranging at least over Z=0.2- 
1 Zq for Lyman Break Galaxies (Kobulnicky & Koo 2000; Pettini et al. 2000, 2001; Teplitz et 
al. 2000), and extending to even lower values for damped Lya systems (e.g., Pettini et al. 1999; 
Prochaska, Gawiser, & Wolfe 2001). To study the effect that metallicity has on the fidelity of 
photometric redshifts, the standard grid of constantly star-forming, Z=0.2 Z Q , dust-free templates 
was used to fit "observed" galaxies that have the same star formation history and dust content as 
the grid, but which have either higher (Zj n = Z®), or lower (Zj„=0.02 Z®) metallicities. 

Fig. 8 shows the results of the fit for "observed" galaxies of solar metallicity (Fig. 9 shows fit 
results for Zj n =0.02 Zq). As Fig. 8 illustrates, photometric redshifts of "observed" galaxies that are 
more metal-rich than the fitting template grid are not strongly affected. At high redshift (zj n ~4) 
the ages derived from the fit are overestimated — as expected given that metal-rich galaxies have 
a stronger 1.6/um bump that makes their age-sensitive A res t<1.6/mi SEDs appear more depressed 
and, hence, older. 

Fit results for "observed" galaxies that are more metal-poor than the template grid are shown 
in Fig. 9. Not surprisingly, at Zi n >3 fit ages are younger than the true ages, agej n . For intrinsically 
old galaxies (agej„=lGyr) this does not affect the recovery of redshift significantly. However, for 
intrinsically young objects (agej n =50Myr), this underestimate of the true age pushes the fit into 
the region at age/ji~10Myr that can be heavily degenerate with z/it<1.5 solutions. Thus, while 
the redshifts of metal-rich high-z galaxies should be recovered faithfully, IRAC colors alone are not 
sufficient to unambiguously identify high-z, very metal-poor galaxies that are dominated by young 
stellar populations (see also Simpson & Eisenhardt, 1999). 

4. SUMMARY AND DISCUSSION 

This paper considers the viability of the 1.6/xm bump as a photometric redshift indicator. In 
§2.2 it was shown that JKHL photometry can be used to estimate the redshifts of z<l galaxies, 
thus demonstrating with real data that the 1.6/xm bump is a viable photometric redshift indicator. 
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The analysis of the dependence of galaxy colors on redshift (§2.3) leads to the conclusion that 
photometric redshift signal at a given redshift comes primarily from two filters: an "anchor" filter 
which at all redshifts of interest remains in the invariant part of the spectrum redward of the 1.6/im 
bump, and a second filter whose redshifted bandpass corresponds to 1.6/um. To cover an extensive 
redshift range and to reduce redshift degeneracies, one needs — in addition to the long- wavelength 
anchor filter that is always redward of (l+z)1.6 / um — multiple filters that span the possible range of 
redshifted 1.6/xm. Filters that are located blueward of the redshifted 1.6//m bump do not constrain 
redshift well, but instead provide information about the age of the stellar population. In this 
manner, the HKL filter set was used to determine redshifts of galaxies at z~0-l (the J filter does 
not provide redshift information but could be used as an age indicator), while the SIRTF IRAC 
3.6/xm, 4.5^m, 5.8/Ltm, 8/Ltm filters can be expected to work at higher redshifts. 

Simulations (§3.1) suggest that SIRTF IRAC colors alone may be sufficient to measure pho- 
tometric redshifts of high-z galaxies, provided that care is taken to reject the degenerate objects 
whose fits claim them to be either at z<1.5 or younger than ~20 Myr (these redshift and age cut-offs 
can be adjusted depending on the levels of completeness and contamination that are acceptable). 
At higher redshifts (z>2.5) low-z contaminants can be separated from young high-z galaxies by 
the addition of UV-optical data that are sensitive to the presence of the Lyman break; at z<1.5, 
the accuracy of photometric redshifts would be restored with the addition of H and K bands. 

As was illustrated in §3.2, the ability to determine photometric redshifts from IRAC photom- 
etry should not be significantly compromised by the unknown metallicities, dust content, and star 
formation histories of high-z galaxies. Provided the caveats about rejecting galaxies with Zfu<1.5 
or agejj t <20Myr are followed, SIRTF-based photometric redshifts should provide a powerful way 
to estimate the redshifts of distant galaxies. With SIRTF IRAC data alone, photometric redshifts 
can be expected to be most accurate for galaxies at z~2, with the accuracy decreasing to higher 
redshifts; they can be expected to be completely inaccurate for galaxies at z<1.5. The epoch at 
zk,2 is the transitional era where the cosmic star formation appears to plateau (e.g., Sawicki, Lin, & 
Yee 1997; Steidel et al. 1999) after its rapid rise with lookback time over 0<z<l (Lilly et al. 1996); 
the expected high accuracy of IRAC photometric redshifts at z^2 means that SIRTF surveys are 
well-poised to study galaxies at this hitherto unexplored epoch of transition. 

Although the 1.6 /jm bump can be used to provide a good estimate of redshift from SIRTF 
data, particularly at zm2, ages of stellar populations cannot be determined with great precision 
from SIRTF data alone at any redshift: large scatter in agefu can be expected in the presence of 
even small photometric uncertainties (§3.1), while significant systematic effects will be present due 
to the unknown metallicity, star formation history, and dust content of the target galaxy (§3.2). 
The addition of optical and near-IR data could alleviate some of these problems: while SIRTF data 
are used to estimate redshifts, near-IR data that span the age-sensitive rest-frame 4000A break 
can be added to estimate ages, and dust effects can be constrained in the rest-UV with optical 
data (see, e.g., Sawicki & Yee, 1998). The addition of spectroscopic redshifts might then free the 
SIRTF photometry to be used to better constrain star formation histories and/or metallicities of 
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distant galaxies. It is beyond the scope of this paper to explore what constraints can be placed 
on the evolutionary parameters of high-z galaxies by multiwavelength photometry augmented by 
spectroscopic redshifts; however, such studies that combine spectroscopy with photometry that 
ranges from rest-frame UV to IR are sure to be an important path to understanding the nature of 
distant galaxies. 

Although §2.2 illustrates that reliable photometric redshifts are possible using the 1.6/im bump, 
the essential test of the 1.6/zm bump as a photometric redshift indicator for SIRTF will come 
only in the form of a comparison between spectroscopic redshifts and photometric ones based on 
SIRTF data. The 1.6^m bump technique is expected to be most useful at z«2 — precisely where 
our ability to determine spectroscopic redshifts is limited by the lack of readily accessible strong 
spectral features. However, we should be able to reliably assemble large samples of z~2 galaxies by 
first using photometric redshifts to select z«2 candidates, then obtaining spectroscopic redshifts 
using deep spectroscopy that targets absorption lines that are weak but accessible in the optical 
(e.g., C IV 1550, Mg II 1670), and finally using the spectroscopic redshifts to test and calibrate the 
photometric redshifts. The 1.6/xm bump photometric redshift technique applied to SIRTF (and, 
later, NGST) data appears to be ideally suited to allow a detailed study of galaxies at the crucial 
epoch around z&2. 

This work has benefitted from my interactions with Doug Johnstone, Keith Matthews, Gerry 
Neugebauer, Tom Soifer, Howard Yee, and the anonymous referee. I thank all of them for their 
diverse contributions. 
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A rest(^ m ) 

Fig. 1. — Model spectral energy distributions from the 1996 version of the GISSEL library of 
Bruzual &; Chariot (1993). Shown are SEDs of solar metallicity stellar populations that are forming 
stars at a constant rate with a Salpeter IMF. The SEDs are normalized at 1.6/um. The 1.6^m bump 
is a prominent feature in all but the very youngest stellar populations. At the bottom of the plot 
are shown the filter transmission curves of the SIRTF IRAC filters (redshifted to show them at 
z=2.7 in the rest-frame of the SEDs) and of the JHKL filter set (here redshifted to z=0A) that 
was used in fitting the nine HDF galaxies in Sec. 2.2. 
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spectroscopic redshift, z spec 

Fig. 2. — Comparison of photometric and spectroscopic redshifts for galaxies in the Hubble Deep 
Field that have JHKL photometry. The photometric redshifts were derived using only the JHKL 
data and so rely solely on the 1.6^m bump. The scatter between photometric and spectroscopic 
redshifts is small (o"Az = 0.09) demonstrating that the 1.6/txm bump is a viable redshift indicator. 
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Fig. 3. — Dependence of infrared galaxy colors 
on redshift. Colors were computed from solar 
metallicity, constant star formation rate GISSEL 
models. Within each panel, the solid lines are, 
from top to bottom, for model ages of 20 Gyr, 
10 Gyr, 2 Gyr, 1 Gyr, 100 Myr, and 20 Myr, 
and the dotted line is for the extremely young 1 
Myr-old model. 
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Fig. 4. — Redshift dependence of galaxy colors in 
the SIRTF IRAC bandpasses. As in Fig. 3, solid 
lines show colors of solar metallicity, constant 
star formation rate GISSEL spectral synthesis 
models spanning ages of 20 Myr - 20 Gyr, while 
the dotted line is for the extremely young, 1 Myr- 
old model. 
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Fig. 5. — Confidence in recovery of redshift 
from SIRTF IRAC photometry. Constantly star- 
forming galaxies of 20% solar metallicity, and 
with ages of 50 Myr (left-hand panels) and 1 Gyr 
(right-hand panels) and located at z=0.5, 1, 2, 
3, 4 (from top to bottom), were fitted with a 
grid of templates with the same metallicity and 
star formation history. Filled circles indicate the 
best fit solutions (zf it , agefu) which are, by con- 
struction, identical to the true, input values of 
redshift and age (zi n , agej n ). The contours en- 
close regions of parameter space that will contain 
95% of objects which have (0.05, 0.05, 0.075, 0.1) 
mag uncertainties in the four IRAC filters. 
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Fig. 6. — The effect of mismatch between the 
amounts of dust in the fitted galaxy and the fit- 
ting grid of templates. Asterisks show the input 
redshifts and ages (zi n , agej n ), while filled circles 
indicate the best fit solutions (zfu, age/j t ). The 
"observed" objects are enshrouded in substan- 
tial amounts of dust (E(B — V)=0.3), while the 
fitting templates are dust-free. The close match 
between Z{ n and zju values in this Figure, as 
well as the similarity in the shape of confidence 
regions between this Figure and Fig. 5, indicates 
that the recovery of redshift information is not 
strongly affected by a mismatch in dust obscu- 
ration between "observed" galaxies and fitting 
templates. 
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Fig. 7. — The effect of mismatch between the 
star formation histories of the fitted galaxy and 
the fitting templates. The meaning of symbols 
and contours is the same as in Figs. 5 & 6. The 
"observed" objects have been formed in instanta- 
neous bursts of star formation, while the fitting 
templates are made from models with constant 
star formation. The assumed star formation his- 
tory does not appear to have a strong effect on 
the ability to recover the correct redshift. 
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Fig. 8. — The effect of mismatch in the metal- 
licities of the object and the fitting grid of tem- 
plates. The meaning of symbols and contours is 
the same as in Figs. 5-7. The "observed" object 
has solar metallicity, while the fitting grid has 
Z=0.2 Z Q . 
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Fig. 9. — The effect of mismatch in the metal- 
licities of the object and the fitting grid of tem- 
plates. The meaning of symbols and contours is 
the same as in Figs. 5-8. The "observed" object 
has a metallicity of 0.02 Z Q , while the fitting 
template grid has Z=0.2 Z@. 
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Table 1. Photometry and redshifts of HDF galaxies detected in L-band 





R.A. a 


Dec. a 


Jab° 


H A B h 


K AB h 


L3.2fmi,AB h,C 


Zspec 


Zphot 


12 


36 


44.017 


+62 


12 


50.11 


20.71+0.09 


20.33+0.06 


20.03+0.06 


20.04+0.14 


0.557 


0.70 


12 


36 


44.626 


+62 


13 


04.29 


20.78+0.09 


20.34+0.06 


20.10+0.06 


20.24+0.13 


0.485 


0.55 


12 


36 


48.088 


+62 


13 


09.21 


20.10+0.07 


19.76+0.05 


19.56+0.05 


19.74+0.09 


0.476 


0.50 


12 


36 


49.435 


+62 


13 


46.92 


18.09+0.03 


17.90+0.02 


17.99+0.02 


18.77+0.04 


0.089 


0.15 


12 


36 


49.541 


+62 


14 


06.85 


21.34+0.12 


20.96+0.09 


20.59+0.07 


20.55+0.17 


0.752 


0.80 


12 


36 


51.783 


+62 


13 


53.85 


20.82+0.10 


20.36+0.06 


20.06+0.06 


20.27+0.14 


0.557 


0.55 


12 


36 


53.916 


+62 


12 


54.26 


20.58+0.09 


20.17+0.06 


19.87+0.05 


20.10+0.10 


0.642 


0.55 


12 


36 


55.460 


+62 


13 


11.63 


21.10+0.11 


20.65+0.07 


20.26+0.06 


20.27+0.17 


0.968 


0.80 


12 


36 


56.675 


+62 


12 


45.51 


20.13+0.07 


19.75+0.05 


19.52+0.05 


19.80+0.08 


0.518 


0.45 



a From Hogg et al. (2000). 

b Magnitudes within 2"diameter circular apertures. 

c Based on "total" magnitudes of Hogg et al. (2000) and converted from Vega to AB normalization 
assuming Ls,2fim,AB=Ls ,Vega+2.77. 



